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We study the effect of investor inertia on stock price fluctuations with a market microstructure model comprising 
many small investors who are inactive most of the time. It turns out that semi-Markov processes are tailor 
made for modelling inert investors. With a suitable scaling, we show that when the price is driven by the 
market imbalance, the log price process is approximated by a process with long range dependence and non- 
Gaussian returns distributions, driven by a fractional Brownian motion. Consequently, investor inertia may lead 
to arbitrage opportunities for sophisticated market participants. The mathematical contributions are a functional 
central limit theorem for stationary semi-Markov processes, and approximation results for stochastic integrals of 
continuous semimartingales with respect to fractional Brownian motion. 
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1. Introduction and Motivation We prove a functional central limit theorem for stationary semi- 
Markov processes in which the limit process is a stochastic integral with respect to fractional Brownian 
motion. Our motivation is to develop a probabilistic framework within which to analyze the aggregate 
effect of investor inertia on asset price dynamics. We show that, in isolation, such infrequent trading 
patterns can lead to long-range dependence in stock prices and arbitrage opportunities for other more 
"sophisticated" traders. 

o 

1.1 Market Microstructure Models for Financial Markets In mathematical finance, the dy- 
namics of asset prices are usually modelled by trajectories of some exogenously specified stochastic process 
defined on some underlying probability space P). Geometric Brownian motion has long become 

the canonical reference model of financial price fluctuations. Since prices are generated by the demand of 
market participants, it is of interest to support such an approach by a microeconomic model of interacting 
agents. 



In recent years there has been increasing interest in agent-based models of financial markets. These 
models are capable of explaining, often through simulations, many facts like the emergence of herding 
behavior [41], volatility clustering [42] or fat-tailed distributions of stock returns [17] that are observed 
in financial data. Brock and Hommes [TUl E] proposed models with many traders where the asset price 
process is described by deterministic dynamical systems. From numerical simulations, they showed that 
financial price fluctuations can exhibit chaotic behavior if the effects of technical trading become too 
strong. 

Follmer and Schweizer 27J took the probabilistic point of view, with asset prices arising from a sequence 
of temporary price equilibria in an exogenous random environment of investor sentiment; see |25) . |32j or 
[21)] for similar approaches. Applying an invariance principle to a sequence of suitably defined discrete 
time models, they derived a diffusion approximation for the logarithmic price process. Duffie and Protter 
[22j also provided a mathematical framework for approximating sequences of stock prices by diffusion 
processes. 

All the aforementioned models assume that the agents trade the asset in each period. At the end of 
each trading interval, the agents update their expectations for the future evolution of the stock price and 
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formulate their excess demand for the following period. However, small investors are not so efficient in 
their investment decisions: they are typically inactive and actually trade only occasionally. This may be 
because they are waiting to accumulate sufficient capital to make further stock purchases; or they tend 
to monitor their portfolios infrequently; or they are simply scared of choosing the wrong investments; or 
they feel that as long-term investors, they can defer action; or they put off the time-consuming research 
necessary to make informed portfolio choices. Long uninterrupted periods of inactivity may be viewed 
as a form of investor inertia. The focus of this paper is the effect of such investor inertia on asset prices 
in a model with asynchronous order arrivals. See [37] for an alternative micro-structure model with 
asynchronous trading. 

1.2 Inertia in Financial Markets Investor inertia is a common experience and is well documented. 
The New York Stock Exchange (NYSE)'s survey of individual shareownership in the United States, 
"Shareownership2000" [46] . demonstrates that many investors have very low levels of trading activity. 
For example they find that "23 percent of stockholders with brokerage accounts report no trading at all, 
while 35 percent report trading only once or twice in the last year" (see pages 58-59). The NYSE survey 
(e.g. Table 28) also reports that the average holding period for stocks is long, for example 2.9 years in 
the early 90 's. 

Empirical evidence of inertia also appears in the economic literature. For example, Madrian and Shea 
[35] looked at the reallocation of assets in employees' individual 401 (k) (retirement) plansQ and found "a 
status quo bias resulting from employee procrastination in making or implementing an optimal savings 
decision." A related study by Hewitt Associates (a management consulting firm) found that in 2001, four 
out of five plan participants did not do any trading in their 401 (k)s. Madrian and Shea explain that "if 
the cost of gathering and evaluating the information needed to make a 401 (k) savings decision exceeds 
the short-run benefit from doing so, individuals will procrastinate." The prediction of Prospect Theory 
[55] that investors tend to hold onto losing stocks too long has also been observed ([50]). 

A number of microeconomic models study investor caution with regard to model risk, which is termed 
uncertainty aversion. Among others, Dow and Werlang ([H]) and Simonsen and Werlang ([H]) considered 
models of portfolio optimization where agents are uncertain about the true probability measure. Their 
investors maximize their utility with respect to nonadditive probability measures. It turns out that 
uncertainty aversion leads to inertia: the agents do not trade the asset unless the price exceeds or falls 
below a certain threshold. 

We provide a mathematical framework for modelling investor inertia in a simple microstructure model 
where asset prices result from the demand of a large number of small investors whose trading behavior 
exhibits inertia. To each agent a, we associate a stationary semi-Markov process x a — (a;")t>o on a 
finite state space which represents the agent's propensity for trading. The processes x a have heavy-tailed 
sojourn times in some designated "inert" state, and relatively thin-tailed sojourn times in various other 
states. Semi-Markov processes are tailor made to model individual traders' inertia as they generalize 
Markov processes by removing the requirement of exponentially distributed, and therefore thin-tailed, 
holding times. In addition, we allow for a market- wide amplitude process ^, that describes the evolution 
of typical trading size in the market. It is large on heavy-trading days and small on light trading days. 
We adopt a non-Walrasian approach to asset pricing and assume that prices move in the direction of 
market imbalance. We show that in a model with many inert investors, long range dependence in the 
price process emerges. 

1.3 Long Range Dependence in Financial Time Series The observation of long range de- 
pendence (sometimes called the Joseph effect) in financial time series motivated the use of fractional 
Brownian motion as a basis for asset pricing models; see, for instance, [44] or [19] , By our invariance 
principle, the drift-adjusted logarithmic price process converges weakly to a stochastic integral with re- 
spect to a fractional Brownian motion with Hurst coefficient H > h. Our approach may thus be viewed 
as a microeconomic foundation for these models. A recent paper that proposes entirely different eco- 
nomic foundations for models based on fractional Brownian motion is [36] . An approximation result for 

1 A 401k retirement plan is a special type of account funded through pre-tax payroll deductions. The funds in the account 
can be invested in a number of different stocks, bonds, mutual funds or other assets, and are not taxed on any capital gains, 
dividends, or interest until they are withdrawn. The retirement savings vehicle was created by United States Congress in 
1981 and gets its name from the section of the Internal Revenue Code that describes it. 
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fractional Brownian motion in the context of a binary market model is given in [52] . 

As is well known, fractional Brownian motion processes are not semimartingales, and so these models 
may theoretically allow arbitrage opportunities. Explicit arbitrage strategies for various models were con- 
structed in [JS] , [T3] and [5] . These strategies capitalize on the smoothness of fractional Brownian motion 
(relative to standard Brownian motion) and involve rapid trading to exploit the fine-scale properties of 
the process' trajectories. As a result, in our microstructure model, arbitrage opportunities may arise for 
other, sufficiently sophisticated, market participants who are able to take advantage of inert investors by 
trading frequently. We discuss a simple combination of both inert and active traders in Section 12.31 

Evidence of long-range dependence in financial data is discussed in |19j . Bayraktar et al. |S] studied 
an asymptotically efficient wavelet-based estimator for the Hurst parameter, and analyzed high frequency 
S&P 500 index data over the span of 11.5 years (1989-2000). It was observed that, although the Hurst 
parameter was significantly above the efficient markets value of H — i up through the mid-1990s, it 
started to fall to that level over the period 1997-2000 (see Figure [1]). They suggested that this behavior 
of the market might be related to the increase in Internet trading, which is documented, for example, in 
NYSE's Stockownership2000 [46], pQ, and [Hj, who find that "after 18 months of access, the Web effect 
is very large: trading frequency doubles." Barber and Odean [2] find that "after going online, investors 
trade more actively, more speculatively and less profitably than before" . Similar empirical findings were 
recently reached, using a completely different statistical technique in [B]. Thus, the dramatic fall in the 
estimated Hurst parameter in the late 1990s can be thought of as a posteriori validation of the link our 
model provides between investor inertia and long-range dependence in stock prices. 
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Figure 1: Estimates of the Hurst exponent of the S&P 500 index over 1990s, taken from Bayraktar, Poor and 
Sircar (\^). 

We note the evidence of long memory in stock price returns is mixed. There are several papers in the 
empirical finance literature providing evidence for the existence of long memory, yet there are several 
other papers that contradict these empirical findings; see e.g. [5] for an exposition of this debate and 
references. However, long memory is a well accepted feature in volatility (squared and absolute returns) 
and trading volume (see e.g. [18] and [20]). The mathematical results of this paper might also be seen 
as an intermediate step towards a microstructural foundation for this phenomenon. 
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1.4 Mathematical Contributions We establish a functional central limit theorem for semi-Markov 
processes (Theorem[2j]below) which extends the results of Taqqu et. al. [54], who proved a result similar 
to ours for on/off processes, that is, semi-Markov processes taking values in the binary state space {0, 1}. 
Their arguments do not carry over to models with more general state spaces. Our approach builds on 
Markov renewal theory. We also demonstrate (see Example 13. ip that there may be a different limit 
behavior when the semi-Markov processes are centered, a situation which cannot arise in the binary 
case. Taqqu and Levy 53J considered renewal reward processes with heavy tailed renewal periods and 
independent and identically distributed rewards. They assume a general state space, but the distributions 
of the length of renewal periods does not depend on the current state; for an extension to the case of 
heavy-tailed rewards, see [40] . A recent paper [45] studies the binary case under a different limit taking 
mechanism; see also [28] . 

Binary state spaces are natural for modelling internet traffic, but for many applications in Economics 
or Queueing Theory, it is clearly desirable to have more flexible results that apply to general semi-Markov 
processes on finite state spaces. In the context of a financial market model, it is natural to allow for both 
positive (buying), negative (selling) and a zero (inactive) state. Our results also have applications to 
complex multi-level queueing networks where the level-dependent holding-time distributions are allowed 
to have slowly decaying tails. They may serve as a mathematical basis for proving heavy-traffic limits in 
the network models studied in, e.g. [53], [23] and [55] . 

We allow for limits which are integrals with respect to fractional Brownian motion proving an approx- 
imation result for stochastic integrals of continuous semimartingales with respect to fractional Brownian 
motion. Specifically, we consider a sequence of good semimartingales {^> n } and a sequence of stochastic 
processes {A^™} having zero quadratic variation and give sufficient conditions which guarantee that joint 
convergence of (A™, \E ,n ) to [B , VP), where B H is a fractional Brownian motion process with Hurst pa- 
rameter H > i, and "J is a continuous semimartingale, implies the convergence of the stochastic integrals 
J \& n dX n to J ^dB H . In addition, we obtain a stability result for the integral of a fractional Brownian 
motion with respect to itself. These results may be viewed as an extension of Theorem 2.2 in [35] beyond 
the semimartingale setting. 

The remainder of this paper is organized as follows. In Section 2, we describe the financial market model 
with inert investors and state the main result. Section 3 proves a central limit theorem for stationary 
semi-Markov processes. Section 4 proves an approximation result for stochastic integrals of continuous 
semimartingales with respect to fractional Brownian motion. 

2. The microeconomic setup and the main results We consider a financial market with a set 
A := {ai,a2, . . . ,&at} of agents trading a single risky asset. Our aim is to analyze the effects investor 
inertia has on the dynamics of stock price processes. For this we choose the simplest possible setup. 
In particular, we model right away the behavior of individual traders rather than characterizing agents' 
investment decisions as solutions to individual utility maximization problems. Such an approach has also 
been taken in [25], [57], gT], [25] and [37] for example. 

We associate to each agent a G A a continuous-time stochastic process x a = (xf) t > on a finite state 
space E, containing zero. This process describes the agent's trading mood. He accumulates the asset at 
a rate VPtx" at time t > 0. The random quantity ty t > describes the size of a typical trade at time t, 
and may be negative, indicating the agent is selling. Agents do not trade at times when xl = 0. We 
therefore call the state the agents' inactive state. 

Remark 2.1 In the simplest setting, x a G { — 1,0,1}, so that each investor is either buying, selling or 
inactive, and V? = 1: there is no external amplification. Even here, the existing results in \54^ do not 
apply because the state space is not binary. 

The holdings of the agent a G A and the "market imbalance" at time t > are given by 



respectively. Hence the process (J t Af ) t >o describes the stochastic evolution of the market imbalance. In 
our microstructure model, market imbalance will be the only component driving the dynamics of asset 
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prices. All the orders are received by a single market maker who clears the trades and sets prices as to 
reflect the incoming order flows. That is, the market maker sets prices in reaction to the evolution of 
market imbalances. 

Remark 2.2 In our continuous time model buyers and sellers arrive at different points in time. Hence 
the economic paradigm that a Walrasian auctioneer can set prices such that the markets clear at the end 
of each trading period does not apply. Rather, temporary imbalances between demand and supply will 
occur, and prices are assumed to reflect the extent of the current market imbalance. In the terminology 
explained in \29^ . ours is a model of a "continuous market (trading asynchronously during continuous 
intervals of time)", rather than a "call market (trading synchronously at pre-established discrete times)". 
As Garman reports, the New York Stock Exchange was a call market until 1811, and since then has 
become a continuous market. (See also Chapter 1 of \41j -) 

We consider the pricing rule 

dSf = J2^t^dt and so Sf = 5 + if , (2) 

aeA 

for the evolution of the logarithmic stock price process S N — (S^)t>o- This is the simplest mechanism 
by which incoming buy orders increase the price and sell orders decrease the price. Other choices might 
be utilized in a future work studying, for example, the effect of a nonlinear market depth function, but 
these are beyond the scope of the present work. (The choice of modelling the log-stock price is simply 
standard finance practice to define a positive price). 

Kruk (|37j) considered a model for continuous auction market, in which order arrivals are modelled by 
independent renewal processes. There are a finite number of possible prices, and agents randomly submit 
price dependent limit orders. These are stored in the order book waiting the arrival of matching orders. 
Kruk finds a limiting distribution of the outstanding number of buy/sell orders at one of the possible 
prices. In contrast, our aim is to find the limiting price process that is driven by the market imbalance 
under different assumptions on the market micro-structure. 

2.1 The dynamics of individual behavior Next, we specify the probabilistic structure of the 
processes x a . We assume that the agents are homogeneous and that all the processes x a and ^ are 
independent. It is therefore enough to specify the dynamics of some reference process x = (xt)t>o- In 
order to incorporate the idea of market inertia as defined by Assumption 12.21 below, we assume that x is 
a semi-Markov process defined on some probability space (f2, J- ', P) with a finite state space E. Here E 
may contain both positive and negative values and we assume G E. The process x is specified in terms 
of random variables : — > E and T n : Q — > R+ which satisfy = To < T± < ■ ■ ■ almost surely and 

P{£n+i = j, T n+ i — T n < t|£i, £ n ; Ti, T n } = P{£ n+1 = j, T n+ i — T n < t\£ n } 

for each n G N, j G E and all t G through the relation 

X t = ^2^nl[T n ,T n+1 )(t). (3) 
n>0 

Remark 2.3 In economic terms, the representative agent's mood in the random time interval [T„,T„ + i) 
is given by £„. The distribution of the length of the interval T n +i — T n may depend on the sequence 
{£n}neN through the states and This allows us to assume different distributions for the lengths 

of the agents' active and inactive periods, and in particular to model inertia as a heavy-tailed sojourn 
time in the zero state. 

Remark 2.4 In the present analysis of investor inertia, we do not allow for feedback effects of prices into 
agents ' investment decisions. While such an assumption might be justified for small, non-professional 
investors, it is clearly desirable to allow active traders' investment decisions to be influenced by asset 
prices. When such feedback effects are allowed, the analysis of the price process is typically confined to 
numerical simulations because such models are difficult to analyze on an analytical level. An exception is 
a recent paper [26] where the impact of contagion effects on the asymptotics of stock prices is analyzed in 
a mathematically rigorous manner. One could also consider the present model as applying to (Internet 
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or new economy) stocks where no accurate information about the actual underlying fundamental value 
is available. In such a situation, price is not always a good indicator of value and is often ignored by 
uninformed small investors. 

We assume that x is temporally homogeneous under the measure P, that is, 

nin+i = j, T n+1 - T n < t\£ n =i} = Q(i,j, t) (4) 

is independent of n G N. By Proposition 1.6 in [TS], this implies that {(„} n6 m is a homogeneous Markov 
chain on E whose transition probability matrix P — (py ) is given by 

Pij = 1™ Q(i,3,t)- 

t— »oo 

Clearly, x is an ordinary temporally homogeneous Markov process if Q takes the form 

Q(i,j,t)= Pij (l-e- Xit ). (5) 
We assume that the embedded Markov chain {£, n }n&i satisfies the following condition. 

Assumption 2.1 For all i,j G E, i ^ j we have that pij > 0. In particular, there exists a unique 
probability measure n on E such that ttP = it . 

The conditional distribution function of the length of the n-th sojourn time, T n+ \ — T n , given £ n +i and 
is specified in terms of the semi-Markov kernel {Q(i, j,t);i, j G E, t > 0} and the transition matrix 
P by 

G(i,j,t) := 91hlA = ¥{T _ Tn < t \ in = t = j} . (6) 

Pij 

For later reference we also introduce the distribution of the first occurrence of state j under P, given 
.To = i. Specifically, for i ^ j, we put 

F(i,j,t):=F{r j <t\x Q = i}, (7) 

where Tj :— ini{t > : xt = j}- We denote by F(j,j,-) the distribution of the time until the next 
entrance into state j and by 

Vj ■= J tF(j,j,dt) (8) 

the expected time between two occurrences of state j E E. Further, we recall that a function L : R + — > M + 
is called slowly varying at infinity if 

L(xt) 

hm — — - = 1 for all x > 

t^oo L(t) 

and that f(t) ~ g(t) for two functions f,g : R+ — » K+ means limt^oo ^jy = 1- 

Assumption 2.2 (i) The average sojourn time at state i G E is finite: 

nii := E[T n+ i - T n \£ n = i] < oo. (9) 

Here E denotes the expectation operator with respect to P. 

(ii) There exists a constant 1 < a < 2 and a locally bounded function L : R + — > R + which is slowly 
varying at infinity such that 

P{T n+1 - T n > i|£„ = 0} - t- a L(t). (10) 
(Hi) The distributions of the sojourn times at state i =/= satisfy 

P{T„+i - T n > t\£ n = *} n 

hm -, — — t— . . = 0. 

f.^co t-i a + l )L(t) 

(iv) The distribution of the sojourn times in the various states have continuous and bounded densities 
with respect to Lebesgue measure on M. + . 
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Our condition (|10p is satisfied if, for instance, the length of the sojourn time at state € E is 
distributed according to a Pareto distribution. Assumption 12.21 (iii) reflects the idea of market inertia: 
the probability of long uninterrupted trading periods is small compared to the probability of an individual 
agent being inactive for a long time. It is stronger than the corresponding assumption for the binary case 
in [54] where the sojourn time in the only other state may in fact be as heavy tailed. For our economic 
application, however, it is natural to think of the sojourn times in the various active states as being thin 
tailed, such as in the exponential distribution, since small investors typically do not trade continually for 
long periods. 

2.2 An invariance principle for semi-Markov processes In this section, we state our main 
results. With our choice of scaling, the logarithmic price process can be approximated in law by the 
stochastic integral of "J with respect to fractional Brownian motion B H where the Hurst coefficient H 
depends on a. The convergence concept we use is weak convergence on the Skorohod space ID of all 
real- valued right continuous processes with left limits. We write C- lim, woo Z n — Z if a sequence of 
D-valued stochastic processes {Z n } n ^i, converges in distribution to Z. 

In order to derive our approximation result, we assume that the semi- Markov process x is stationary. 
Under Assumption ! 2 . 1 1 stationarity can be achieved by a suitable specification of the common distribution 
of the initial state £o and the initial sojourn time T±. We denote the distribution of the stationary semi- 
Markov processes by V* . The proof follows from Theorem 4.2.5 in [SJ, for example. 

Lemma 2.1 In the stationary setting, that is, under the law P* the following holds: 

(i) The joint distribution of the initial state and the initial sojourn time takes the form 

poo 

F>*{£ = k,T x > t} = = - / h{k,s)ds. (11) 



SjG_E n 3""3 Jt 

Here rm denotes the mean sojourn time in state i € E as defined by (OJ), and for i G E, 

h{i,t) = l-Y,Q(h3>*) (12) 

is the probability that the sojourn time at state i G E is greater than t. 
(ii) The law v = (fk)keE of xt in the stationary regime is given by 

Vk = ^ • (13) 

(iii) The conditional joint distribution of (£\,T\), given £o is 

r{fi=j ) r 1 <t|e = *} = -^ [ [i-G(kj,s)}ds. (u) 

™k,j Jo 

Here mtj '■= f£°[l ~ G(k,j, s)]ds denotes the conditional expected sojourn time at state k, given 
the next state is j, and the functions G(k,j, ■) are defined in {(|). 

Let us now introduce a dimensionless parameter e > 0, and consider the rescaled processes x^ £ . For 
e small, x^, e is a "speeded-up" semi-Markov process. In other words, the investors' individual trading 
dispensations are evolving on a faster scale than Observe, however, that we are not altering the main 
qualitative feature of the model. That is, agents still remain in the inactive state for relatively much 
longer times than in an active state. 

Mathematically, there is no reason to restrict ourselves to the case where *S> is non-negative. Hence 
we shall from now on only assume that ^ is a continuous semimartingale. Given the processes & and x a 
(a € {ai, . . . , ajv}), the aggregate order rate at time t is given by 

Yt" = Y,*t X t/e- (15) 

Let n := E*x t and X e ' N = (X^ N ) <t<T (T > 0) be the centered process defined by 

X t' N : = / E ^ ds = I Y s' N ds - f ^ s ds. (16) 

Jo „ cA Jo Jo 
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We are now ready to state our main result. Its proof will be carried out in Sections [3] and [¥] The definition 
of the stochastic integral with respect to fractional Brownian motion will also be given in Section [4j 

Theorem 2.1 Let ^ = (^t)t>o be a continuous semimartingale on (f2, .F, P*) with a decomposition 
= M + A, in which M is a local martingale and A is an adapted process of finite variation. We assume 
that K{[M , M]t} < oo and thatM{\A\x} < oo, where (|^4|*)t>o is the total variation of A. If Assumptions 
7] and \2.2\ are satisfied, and if [iJ^keE > ®> then there exists c > such that the process X e ' N 



satisfies 



£-lim£- lim Xf' N = [c ^sdB") . (17) 



0<t<T 



/ 0<t<T 



Here the Hurst coefficient of the fractional Brownian motion process B H is H = > i . 

Observe that Theorem 12.11 does not apply to the case /i = 0. For centered semi-Markov processes 
x a , Example 13.11 below illustrates that the limiting process depends on the tail structure of the waiting 
time distribution in the various active states. This phenomenon does not arise in the case of binary state 
spaces. 

Remark 2.5 (i) Theorem \2.1\ says the drift- adjusted logarithmic price process in our model of inert 
investors can be approximated in law by the stochastic integral of W with respect to a fractional 
Brownian motion process with Hurst coefficient H > | . 

(ii) In a situation where the processes x a are independent, stationary and ergodic Markov processes 
on E, that is, in cases where the semi-Markov kernel takes the form £5J), it is easy to show that 

£-lim£- lim ( —^X^' N \ = (c [ ^ s dW. 



/ 0<t<T \ JO / 0<t<T 

where (Wt)t>o is a standard Wiener process. Thus, if the market participants are not inert, that 
is, if the distribution of the lengths of the agents' inactivity periods is thin-tailed, no arbitrage 
opportunities emerge because the limit process is a semimartingale. 



The proof of Theorem 12. II will be carried out in two steps. In Section [3] we prove a functional central 
limit theorem for stationary semi-Markov processes on finite state spaces. In Section [4] we combine our 
central limit theorem for semi-Markov processes with extensions of arguments given in [38 to obtain (|17[) . 

2.3 Markets with both Active and Inert Investors It is simple to extend the previous analysis 
to incorporate both active and inert investors. Let p be the ratio of active to inert investors. We associate 
to each active trader b £ {1,2,..., pN} a stationary Markov chain y b = (j/j )t>o on the state space E. The 
processes y b are independent and identically distributed and independent of the processes x a . The thin- 
tailed sojourn time in the zero state of y b reflects the idea that, as opposed to inert investors, these agents 
frequently trade the stock. We assume for simplicity that \& = 1. With Y t e ' N = Ylb=i f^t/e ~ ^*2/o) an d 
Xt' N := f Q Y^' N ds, it is straightforward to prove the following modification of Theorem 12.11 

Proposition 2.1 Let x a (a — 1, 2, ... . N) be semi-Markov processes that satisfy the assumption of The- 
orem [2~l[ If y b (b = 1,2, . . . , pN) are independent stationary Markov processes on E, then there exist 
constants c\ , C2 > such that 

£-lim£- lim ( } Xf' N + -±=X?' N ) = ( Cl Bf + c 2 Jp Wt)^^ . 

eio iv^oo yei-By/NHe-i) 1 ^/We 4 / o<t<T 

Here, W = {Wt)t>o is a standard Wiener process. 



Thus, in a financial market with both active and inert investors, the dynamics of the asset price 
process can be approximated in law by a stochastic integral with respect to a superposition, B H + SW, 
of a fractional and a regular Brownian motion. It is known ([12]) that B H + 8W is a semimartingale for 
any S ^ 0, if H > |, that is, if a < |, but not if H 6 (5, |]. Thus, no arbitrage opportunities arise 



Bayraktar, Horst, Sircar: Investor Inertia 

Mathematics of Operations Research xx(x), pp. xxx-xxx, ©200x INFORMS 



9 



if the small investors are "sufficiently inert." The parameter a can also be viewed as a measure for the 
fraction of small investors that are active at any point in time. Hence, independent of the actual trading 
volume, the market is arbitrage free in periods where the fraction of inert investors who are active on the 
financial market is small enough. 



3. A limit theorem for semi-Markov processes This section establishes Theorem 12.11 for the 

special case *S? = 1. We approach the general case where "J is a continuous semimartingale in Section [U 
Here we consider the situation where 

Y t e ' N = V x a tle and where Xp N = [ Y*> N ds - N/J,t, 

and prove a functional central limit theorem for stationary semi-Markov processes. Our Theorem 13.11 
below extends the results in [54] to situations where the semi-Markov process takes values in an arbitrary 
finite state space. The arguments given there are based on results from ordinary renewal theory, and 
do not carry over to models with more general state spaces. The proof of the following theorem will be 
carried out through a series of lemmas. 

Theorem 3.1 Let H = Under the assumptions of Theorem [Ol 

£-lim£- Hm f , „ ,L, „ X? N ) = (cflf ) 0<{<T . (18) 



Q<t<T 



Let 7 be the covariance function of the semi-Markov process (x t ) t >o under P*, and consider the case 
e = 1. By the Central Limit Theorem, and because x is stationary, the process Y — (Y t )t>o defined by 

Y t = C- lim -^={Yt ,N - iV>) (19) 



is a stationary zero-mean Gaussian process. It is easily checked that the covariance function of the process 
(-^=Yt' N ) is also 7 for any AT, and hence for Yt. By standard calculations, the variance of the aggregate 

process (J* Y s ds) at time t > is given by 

Var(t) :=Var^y Y s dsj=2^ ^ j(u)duj dv. (20) 

In the first step towards the proof of Theorem 13. li we can proceed by analogy with [54] . We are 
interested in the asymptotics as e J, of the process 



XI := / Y s/£ ds, (21) 
Jo 

which can be written X^ = e J^ e Y s ds. Therefore the object of interest is the large t behavior of Var(£). 
Suppose that we can show 

Var(i) - c 2 t 2H L(t) as t -> oo. (22) 
Then the mean-zero Gaussian processes X € = (Xf) t >Q have stationary increments and satisfy 

limE* | 1 Yf I =c 2 t 2H . (23) 

40 \ e i-By/L(FI) V 

Since the variance characterizes the finite dimensional distributions of a mean-zero Gaussian pro- 
cess with stationary increments, we see that the finite dimensional distributions of the process 



— - -, = XI J converge to (cB^ t>Q whenever (|22p holds. The following lemma gives a suffi- 
cient condition for (|22p in terms of the covariance function 7. 



Lemma 3.1 For [2ty) to hold, it suffices that 

7(t) - c 2 H(2H - l)t 2H - 2 L(t) as t -> 00. (24) 



10 
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Proof. By Proposition 1.5.8 in [5], every slowly varying function L which is locally bounded on 
satisfies 

if /3 > — 1. Applying this proposition to the slowly varying function 

lit) 



we conclude 



' 2 

j(u)dudv ~ y< 2ff L(i), 



10 JO 

and so our assertion follows from (|2"0)) . □ 

Before we proceed with the proof of our main result, let us briefly consider the case [i = which is 
not covered by our theorem. For semi-Markov processes whose "heavy-tailed state" happens to be the 
mean, the structure of the limit process depends on the distribution of the sojourn times in the various 
active statesH 

Example 3.1 We consider the case E = { — 1, 0, 1}, and assume that p-ifi = Pxfi — 1 an d that po.-i — 
Po,i — \ - With v\ = ¥*{x t = 1} > 0, we obtain 

7(£) = v\ {K*\x t xo\xo = 1] + E*[x t x \xo = -1]) • 

Suppose that the sojourns in the inactive state are heavy tailed, and that the waiting times in the active 
states are exponentially distributed with parameter 1. In such a symmetric situation 

E*[x t x \x = ±1] = P*{Ti > t\x Q = ±1} = e"*. 

Therefore, "/(t) = 2v\e~ l . In view of i23\), this yields c > such that 

C-lhnC- lim (^=XZ< N ) =(cW t ) Q<t<T 

elO N— too \ W£j\J 



0<t<T 



for some standard Wiener process W . 



In order to prove Theorem 13. li we need to establish ([24]) . For this, the following representation of 
the covariance function turns out to be useful: in terms of the marginal distribution vi = P*{xt = i} 
(i G E) of the stationary semi-Markov process given in Lemma 12.11 (i) , and in terms of the conditional 
probabilities 

P:(i,j):=V*{x t =j\x =i}, 

we have 

7 (t)= Yl im(Pnhj)-^)- (25) 

It follows from Proposition 6.12 in [15] , for example, that P t * (i, j) — * Vj as t — > co. Hence limt^oo j(t) = 0. 
In order to prove Theorem 13. li however, we also need to show that this convergence is sufficiently slow. 
We shall see that the agents' inertia accounts for the slow decay of correlations. It is thus the agents' 
inactivity that is responsible for that fact that the logarithmic price process is not approximated by 
a stochastic integral with respect to a Wiener process, but by an integral with respect to fractional 
Brownian motion. 

We are now going to determine the rate of convergence of the covariance function to 0. To this end, we 
show that P*(i,j) can be written as a convolution of a renewal function with a slowly decaying function 
plus a term which has asymptotically, i.e., for t — > oo, a vanishing effect compared to the first term; 
see Lemma 13.21 below. We will then apply results from [31] and [34] to analyze the tail structure of the 
convolution term. 

Let 

R(i,j, t) := E I % n =j,T„<t} | xa = $a=i\ 
U=o J 

2 We thank Chris Rogers for Example 13. II 
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be the expected number of visits of the process (xt)t>o to state j up to time t in the non-stationary 
situation, i.e., under the measure P, given xq — i. For fixed i,j € E, the function t \— ► R(i,j,t) is a 
renewal function. If, under P, the initial state is j, then the entrances to j form an ordinary renewal 
process and 

oo 

R(j,j ) t) = Y,F n (jJ ) t). (26) 

71 = 

Here F(j,j, •) denotes the distribution of the travel time between to occurrences of state j £ E as defined 
in 0, and F n (j,j, •) is the n-fold convolution of F(j,j, •). On the other hand, if i ^ j, the time until 
the first visit to j has distribution F(i,j, •) under P which might be different from F(j,j, ■). In this case 
R(i,j, •) satisfies a delayed renewal equation, and we have 

R(i,j,t) = f R(j,j,t-u)F(i,j,du). (27) 
Jo 

We refer the interested reader to [15] for a survey on Markov renewal theory. 

Let us now return to the stationary setting and derive a representation for the expected number 
R*(i,j, t) of visits of the process {xt)t>o to state j up to time t under P*, given xq = i. To this end, we 
denote by F*(i,j, •) the distribution function in the stationary setting of the first occurrence of j, given 
xq = i and put 

P t (i,j) :=F{x t =j\x =i}. 
Given the first jump time T\ and given that xt x = i we have that 

F*{x t =j\x Tl =i} = P t -T 1 (i,j) on {t > Ti}. (28) 

Thus, 

R*(i,j,t)= [ R(j,j,t-u)F*(i,j,du). (29) 
Jq 

3.1 A representation for the conditional transition probabilities In this section we derive a 
representation for P t * (i, j) which will allow us to analyze the asymptotic behavior of P t * — Vj . To this 
end, we recall the definition of the joint distribution of the initial state and the initial sojourn time and 
the definition of the conditional joint distribution of given £o from (TTTj) and (TH| respectively. 

We define 

s(i,t) :=P*{£o =i,T 1 >t} and s(i, j, t) = P*{£i = j,Tx < i|£ = *}. (30) 
In terms of these quantities, the transition probability P£(i,j) can be written as 

= + £ [ P t -u(k,j)s(i,k,du). (31) 

Here the first term on the right-hand-side of (|3"Tj) accounts for the P*-probability that Xq = i and that 
the state i survives until time t. The quantity Pt- U (k,j)s(i, k, du) captures the conditional probability 
that the first transition happens to be to state k before time t, given £o = i- Observe that we integrate 
the conditional probability P f _ u (i, j) and not P*_ u (i,j): conditioned on the value of semi-Markov process 
at the first renewal instance the distributions of (x t )t>o under P and P* are the same; see (|2"5|) . 

In the sequel it will be convenient to have the following convolution operation: let h be a locally 
bounded function, and F be a distribution function both of which are defined on R + . The convolution 
F * h of F and h is given by 

F * h(t) := [ h(t - x)F(dx) for * > 0. (32) 
Jo 

Remark 3.1 Since F*h is locally bounded, the map t <—> G * (F * h)(t) is well defined for any distribution 
G on R + . Moreover, G * (F * h)(t) = F*(G*h)(t) = {G * F)*h{t). In this sense distributions acting on 
the locally bounded function can commute. Thus, for the renewal function R = 'Y^ =0 F n associated to F , 
as defined in i26]) , the integral R* h(t) is well defined and R* (G * h)(t) = G* (R* h) (t) = (R*G) * h(t). 



12 
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We are now going to establish an alternative representation for the conditional probability P t *(i, j) 
that turns out to be more appropriate for our subsequent analysis. 

Lemma 3.2 In terms of the quantities s(i,t) and h(i,t) in ilty) and R*(i,j,t), we have 

Pt*(hj) = —Sij + [ h(j,t-s)R*(i,j,ds). (33) 



Proof. In view of (|3Tj) . it is enough to show 

R*(i,j,t) * h(j,t) = / Pt-u{k,j)s(i,k,du). 
keE Jo 

To this end, observe first that F*(i,j, t) can be decomposed as 



F*(i,j,t) = s(i,j,t) + V / F(k,j,t-u)s(i,k,du). 



(34) 



Indeed, s(i,j, t) is the probability that the first transition takes place before time t and happens to be to 
state j € E, and 

ft 

F(k,j, t — u)s(i, k, du) = ¥* {x v = j for some v < t, xt x = k\x$ = i}. 

In view of (23]) and ((Ml), Remark lO yields 

R* (i, j, t) * h(j, t) = R(j, J, t) * F* {%, j, t) * h(j, t) 

= R(j, j, t) * s(i, j, t) * h(j, t) + J2 F ( k , h t) * s(i, k, t) * R(j, j, t) * h(j, t). 

Now recall from Proposition 6.3 in [15], for example, that 

p t(hj)= h(j,t- s)R(i,j,ds). 
Jo 

Thus, by also using (f2"T|) we obtain 

R* t) * h(j, t) = s{i, j, t) * P t {j, j) + Rti, 3, t) * F ( k , i> *) * *(*. k ^ t) * Ki, t) 

= s(i, j, t) * P t (j, j) + R ( k > 0, <) * k, t) * h(j, t) 
= Y,s{i,k,t)*P t {k,j). 

k 

This proves our assertion. □ 

3.2 The rate of convergence to equilibrium Now, our goal is to derive the rates of convergence 
of the mappings t i— > s(i,t) to and t i— > R*(i,j,t) * h(j,t) to Vj, respectively. Due to (|2"5]l it is enough 
to analyze the case i,j ^ 0. To this end we shall first study the asymptotic behavior of the map 
t i— > R*(i,j,t). Since R(j,j,-) is a renewal function, we see from (f29|) and (f33|) that the asymptotics of 
Pt(i,j) can be derived as an application of Theorem I A. II essentially if we can show that 

F*(i,j,t)*h(j,t) = o(F(j,j,t)). (35) 

3.2.1 The tail structure of the travel times Let us first deal with the issue of finding the 
convergence rate of F(j,j,t) = 1 — F(j,j,t) to 0. To this end, we introduce the family of random 
variables 

& = {(ef ij ),i,jeE,£ = 0,1,2,...}, (36) 

such that any two random variables in are independent, and for fixed pair (i, j) the random variables 
B\ j have G(i,j,-) as their common distribution function. To ease the notational complexity we assume 
that the law of T n+ i — T n only depends on £ n . We shall therefore drop the second sub-index from the 
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elements of 0. The random variables (9j)e^ are independent copies of the sojourn time in state i. 
We shall prove Lemma 13.41 below under this additional assumption. The general case where T n+ i — T n 
depends both on £ n and £ ra +i can be analyzed by similar means. 

Let Nl' 3 denote the number of times the embedded Markov chain {£„}„gN visits state k € E before it 
visits state j, given £o = *• By definition, Nj' J = with probability one. We denote by N l,J the vector 
of length \E\ with entries Nl' 3 , and by n = {nk)k£E an element of N' E '. Then we have 

[ k^j 1=1 J 

With G(k,t) := P{T„+i - T n < i|£ n = k}, we can rewrite {37]) as 



F(i,j, t ) = E p i^°+EE^> i 



N itj = n } F{N i>j = n} 



(38) 



1 - G(i, t)*J2 * G n " (k, t)P{N lJ = n}. 
n k ^ 3 



Our goal is now to show that 



for j ^ 0. Here, Si.o = if i ^ and 6i t o — 1 otherwise. The first term on the right-hand-side, 
^ n>0 riP{A f Q J = n}, is the expected number of occurrences of state under P before the first visit to 
state j, given £ = x — i. This quantity is positive, due to Assumption 12.11 In order to prove (|3"5|) , we 
need the following results which appear as Lemma 10 in [34] . 

Lemma 3.3 Let F\, F m be probability distribution junctions such that, for all j ^ i, we have Fj(t) — 
o(Fi(t)). Then for any positive integers m, . . . , n m , 

Moreover, for each u > 0, there exists some K u < oo such that 

1-F?(t) ^ K ^ 1 + u ) 

for all t > 0. 

We are now ready to prove (f39|) . 
Lemma 3.4 Under the assumptions of Theorem \3.1\ we have, for j ^ 0. 

v ' n>0 

Proof. Let us first prove that the expected number of occurrences of state before the first return 
to state j occurs is finite. To this end, we put p = min-jjjjj : i,j 6 E,i ^ j} > 0. Since 

P{A/j j = n} < P{£ m £ j for all m < n} < (1 - p) n , 

we obtain 

J2 n ^{K j = n} < X] n (! ~P) n < 00 • ( 4 °) 

n>0 n>0 

Now, we define a probability measure p, on N^ E \ by 

fl{n} = P{AT lJ = n} 

and put 

A n (t) = G(i,t) * G" fc (fc,t). 
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Since — >■ 1 as t — * oc, the first part of Lemma [3^31 yields 

,. 1-4,(0 / 1-G(i,t) t . G~(t,«) \ 

lim r^r~ = lim ^77 ; — — 1M) + Oi n. 

t-00 t^oo I 1-G(0,<) J V t- a L(t) 1 

From the definition of the measure p,, we obtain 

1 - En A n (t)F{N^ = n} _ _ ( 1 - A n (t) 



t- a L(t) Efl \ t~ a L(t) 

and so our assertion follows from the dominated convergence theorem if we can show that 

1-A n (t) rl ,_. 

s T^WT eL{fi) - (41) 



To verify (|4Tj) . we will use the second part of Lemma 13731 For each u > there exists a constant K u such 
that 

Since 

lim = 1( 

t^co t~ a L{t) 

and because we are only interested in the asymptotic behavior of the function t >— > f-l'ify , we may with 
no loss of generality assume that 

1-G(0,*) 

SUP — T" = 1. 

This yields 

T 1 -r^ l - K41+uyi0+5tX '- (42) 

From (gDJ and (J42J) we get 

% {sup < K u (l + «)*•■ f> - P ) fc (l 

Choosing u < we obtain /3 := (1 — p)(l + u) < 1 and so the assertion follows from 



< CO. 

n>0 



□ 



So far, we have shown that F(i,j,t) ~ t a £(i) (Xm>0 T ^{-^o'"' = "1 + ^i,o J for j 7^ 0. In view of 
Lemmas 13.31 and \3A\ the representation (|34l of F*(i,j,t) yields a similar result for the stationary setting. 



Corollary 3.1 For all i,j ^ we Ziave 

lim — 

t^oo t~ a L(t) 



l im = y n P{iV^ = n}. 



n>0 

Proof. Due to (JH3J), (J30l) and JT4J) we can write 



F*(i,j,t) = ^- [ (l-G(i,j,s))ds + J2^- [ F(k,j,t-u)(l-G(i,k,u))du, (43) 



in. 



'<:l 

and therefore 



f (t,j,t) 



jf (l-G(iJ,s))ds + ^p«Yl- ^— J F(k,j,t-u){l-G{i,k,u))du) . (44) 

3 
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We will now show that 

f°°(l-G(i,j,s))ds 

lim ^ ; = if Pij > and i ^ 0. 45 

t^co t- a L(t) y ' 

To this end, we first apply Proposition 1.5.10 in [5]: if g is a function on M + that satisfies g(t) ~ t~@L(t) 
for j3 > 1, then 

^-/3+l 

5 ( s )ds ~ -5-^r£(t). 



/8-r 

Together with Assumption 12.21 (iii) this proposition implies that 

where h(i, •) is the tail of the distribution of the sojourn time in state i, defined in (JX3J) . The representation 

h(i,t) = 1 - X>iG(M,i) = ^Pij(l - G(M,i)), (47) 

along with (|4"5)) implies (|4"S"]) . 

In order to find the decay rate of the remaining terms of (|4"4")h recall first that N 1 ^ 3 is the number of 
visits to state k before reaching state j, given £ = i, and not counting the first one in the case k = i. 

P{N^ =n}= PikHNo' j =n}+ Pi0 F{N l '° = n - 1}, 

and so 

53 n¥{N^ =n} = J2 Pik £ » P {< J = »} + P* E = ( 48 ) 

n>0 fc^j n>0 n>0 

Note that Assumption 12.11 implies Y^ n >o f{^o J = n} = 1. Since J (1 — G(i,j, s))ds is a distribution 
function whose tail is J. (1 — G(i,j,s))ds, Lemmas 13.31 13.41 together with equations (143)) and (|4"8j) 
imply that 

T.k^P^kil- ^^F{k,j,t-u)(l-G{i,k,u))du\ 
t lim : — — = }^ Plk }^nP{N ' J =n}+ PtQ 

^ ' k^j n>0 

= E nP{N^ = n}. 

n>0 

This completes the proof. 

□ 

The next result shows that the first term on the right-hand-side of (|3"3"|) converges to zero sufficiently 
fast. 



Corollary 3.2 For all i ^ we have 



lim \_' ' , = 0. 

t^oo t- a+1 L(t) 



Proof. The proof is an immediate consequence of Corollary 13. II since F*(i,j,t) > s(i,i)/Vj. □ 

3.2.2 The tail structure of R* * h So far, we have analyzed the tail structure of the distribution 
of the travel time between states i and j ^ 0) in the stationary regime. We are now going to study 
the tail structure of R*(i,j,t) * h(j,t). 

Lemma 3.5 There exists Cj > such that, for j ^= 0. 

R*{i,j,t)*h{j,t)-Vj = Cj 
t^S> t- a+1 L(t) a - I 

for alii e E, iyt 0. 



16 



Bayraktar, Horst, Sircar: Investor Inertia 
Mathematics of Operations Research xx(x), pp. xxx-xxx, ©200x INFORMS 



Proof. Let us fix i,j g E, i,j ^ 0. Using the representation (|2"5)) for the function R*(i,j, •), we need 
to show that 

Urn F * ( z ' 3> *) * ^ *) _ ^ - C J 

i — *oo 

It follows from Fubini's theorem that 

F*(i,j,t)*h(j,t)dt 



t~ a+1 L(t) 


a — 1 


p OO pOC 

= / / h(j,t- 


s)l{ s <t}F*(i,j, ds)dt 


Jo Jo 


pOQ pOO 

= 1 1, *<*'- 


s)dt F* (i, j, ds) 






poo 

= / h(j,t)dt 




Jo 




= m h 





where rrij is the mean sojourn time at state j £ E as defined in ©. Suppose now that we can show that 
the continuous non-negative function z defined by 

z(i,j,t) :=F*(i,j,t)*h(j,t) (49) 

satisfies z(i,j,t) = o(F(i,j,t)), is directly Riemann integrableU and of bounded variation, and its total 
variation function z*(i,j, •), defined in ([57)) below, satisfies z*(i,j,t) — 0(t~ a+1 L(t)), where we define 

L{t) = L(t) nP{N^ = n}. 

n>0 

Then Theorem IA. II yields 

^- /* z(i,j,t-s)R(j,j,ds) ~ — mj 2 r a+1 L(t), (50) 
Vj Jo (" - l)Vj 

because R(j, j, t) — J2 n >o ^™ (j> ^ s a renewa l function, F(j, j, •) is nonsingular and F(j, j, t) ~ t~ a L(t) 
for j 7^ 0. By Propositions 5.5 and 6.12 in [IS], Uj = and so (fSTJ)) would yield 

flb-,j,t),F'(Utyqt)- yj . = m, _ = 
i^co r^'i t a - 1)77? ^ 1 ; 

v y ' v n>0 

and hence prove our assertion with 

Cj : ^^ ; , : : { .\;- „}. (51) 

"i n>0 

We will now establish that (i) z(i,j,t) = o(F(i,j,t)), that (ii) z*(i,j,t) = 0(t~ a+1 L(t)) and z(i,j,-) is 
of bounded variation, and that (in) z(i,j, •) is directly-Riemann integrable. 

(i) Since i, j ^ and because p^o > the probability that the semi-Markov process x visits the 
state before it reaches the state j is positive. Thus, it follows from Assumption 12.21 (iii) and 
from Corollary 13.11 that 

h(j,t)=o(F*(i,j,t)). (52) 
Now it follows from (f52"|) and Lemma |3"U1 that 

F* (i, j, t) + z(i, j, t) = l-F* t) * (1 - h(j, t)) ~ F* (i, j, t). (53) 

Hence z(i,j,t) — o(F* (i, j,t)), and so z(i,j,t) — o(F(i,j,t)) by Lemma and Corollarv l3.H 

(ii) Assumption 12.21 (iv) along with (|3"5)) . the representation (|4*3")l of F*(i,j, •) and (|37|) guarantees 
that F*(i,j, •) and /i(j, •) have a bounded continuous densities f*(i,j, ■) and /i'(j, ■), respectively. 
As a result, z(i,j, •) is absolutely continuous with density 

z'(i,j,t) = f*(i,j,t)+ [ h'(j,t-s)r(i,j,s)ds. (54) 



J See Appendix [B] for the definition of directly Riemann integrability. 
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Since h'(j,t- s) < 0, 



\z'(ij,t)\<f*(ij,t)- ti(j,t-s)f*(ij,s)ds. (55) 



In terms of the distribution function g(j, t) := 1 — h(j,t) we have 

dtJo 



h'(j,t-s)f*(i,j,s)ds = ~ f h(j,t- S )f*(i,j,s)ds + r(iJ,t) 

' (56) 
= -Q t (g(j,t)*F*(i,j,t)). 

As a result the total variation function z*(i^j, •) satisfies 

z*(i,j,t)= / \z'{i,j,s)\ds<F*(i,j,t) + [l-g(j,t)*F*{i,j,t)]. (57) 



Thus, z*(i,j,t) < 2 and z(i,j, •) is of bounded variation. Lemma |3~31 and ([52")) together with ([57)) 
give z*(i,j,t) = 0(F*(i,j,t)), and so Corollary O yields z*(i,j,t) = 0(i~ Q+1 £(*)). 

(iii) In order to prove that z is directly Riemann integrable note first that since h > is decreasing 
and by Assumption 12.21 (i) it is integrable, therefore by Lemma IB. II (i) it is directly Riemann 
integrable. Now Lemma IB. II fii) implies that z is directly Riemann integrable. 

□ 

3.2.3 Proof of the central limit theorem for semi-Markov processes We are now ready to 
prove the main result of this section. 



Proof of Theorem 13.11 By we have the representation 



P t *(i,j) = S -^-8 ij + f h(j,t~ s )R*( l ,j,ds) 



Vi JO 

for the conditional probability that xt = j, given xq = i. Due to Corollary 13. 21 and Lemma 13.51 



P?{i,j)-Vj _ Cj 



t^oo t~ a+1 L(t) a~l 



(»,i^o), 



where Cj is given by dHTJ) . With H = ^ it follows from (gSJ) that 

&w^ = -^hs)^ VACi ' (58) 

By Lemma 13.11 this proves the existence of a constant c such that the finite dimensional distributions 
of the processes — - — Xf I converge weakly to the finite dimensional distributions of the 

V el H V L ( £_1 ) / <t<oo 

fractional Brownian motion process cB H as e | 0, and c is given by 

c 2 = ; l — V ijvj^- Y nf{M> j =n} = , l — -uY 2% V n¥{Nfr j = 

2H 1 - H)(2H- 1 ^ ^ 2ff 1 - H)(2H - 1 W ^ J n 2 ^ X 

v n 1 i,jeE '1 n>0 y Jy ' j£E 3 „>o 

In order to establish tightness, we proceed in two steps, 
(i) We first establish the existence of a constant C < oo such that 

Var(t) < Ct 2H L(t) for all t > 0. (59) 
In view of ([22]) . we can choose a sufficiently large T £ N that satisfies 

Var(t) < 2c 2 t 2H L{t) for all t > T. 
In terms of the random variables Y l :~ //*_Twt Y s ds (i = 1,2, ... ,T) we have 



Var ( / Y s ds 



/ T r t/T \ T 

= Var Y / Y s ds = Yl Gov(Y\ Y 1 ). 

\i=l J<i-l)t/T J lJ=1 
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Since the mean zero Gaussian process (J* Y s ds) has stationary increments, the sequence 
Y 1 , Y 2 , . . . , Y T is stationary. Thus, the Holder inequality yields 



Cov(Y\Y ] ) < A/Vax(y i ) v / Vax(FJ) = Var(Y°), 

and so 

Var(t) < T 2 Var for < t < T. 

In view of ([25]) . the function 7 is bounded: H7H00 < 00. Hence for s € [0, 1] the representation 
(j2"0|) shows that 



Var(s) = 2 f f 1 {u)dudv < s^MU < S 2H \h\\ 
Jo Jo 



x because H € [ -, 1 



This yields ([59]) with C := max{2c 2 , T 2_2H ||7||oo} (after putting L = 1 on [0,T]). 
(ii) Let us now denote by Z e = (Zf) the mean zero Gaussian process with stationary increments 
defined by 

Z e t := 1 Xf (0<t<T). 

Due to Theorem 12.3 in [7], the family of stochastic processes (Z £ ) is tight if the following moment 
condition is satisfied for some constants 8 > 1 and C < 00 and for all sufficiently small e: 

E[Zf 2 - Z e ti ] 2 < C\t 2 - h\ 5 for all < h < t 2 < T. (60) 

In view of ([2"T]). we have 

and so it follows from step (i) that 

nzi 2 zif < C (t 2 - tl r L ^^- h) \ (6i) 

Since L is slowly varying, ^fepr^ tends to 1 as e — * 0, and this convergence is uniform in u over 
compact sets ([S], Theorem 1.2.1). Thus, there exists e* > such that 

L( f y < 2 for all < u < T 
if e < £*, and so (JOT)) follows from ([HI]) with C := 2C. □ 



4. Approximating Integrals with respect to fractional Brownian motion In this section 
we prove an approximation result for stochastic integrals which contains Theorem 12.11 as a special case. 
More precisely we give conditions which guarantee that for a sequence of processes {(W™, Z n )}„ 6 N the 
convergence C- lim„_ >00 (*", Z n ) = {^,B H ) implies the convergence L- lim^oo (*™, Z n , J ^ n dZ n ) = 
(^,B H ,J ^dB H ^. We follow the notational convention in [35] that J XdY (without displaying the 
running variable) denotes J X s _dY s . 

All stochastic integrals in this section are understood as limits in probability of Stieltjes-type sums 
which can be described as follows: given a probability space (fi,.F,P) and a filtration {Tt : t > 0} of 
sub-sigma- fields of T consider two adapted stochastic processes <j> and Z; we say that the integral J <fidZ 
exists if for any T < 00 and for each sequence of partitions {t'}; £ n, t 1 = (r\, r|, . . . , tL), of the interval 
[0, T] that satisfies limi^oomaxi |r/ +1 — t\\ = 0, 

/ <p s -dZ s =P- lim y^<p T i(Z T , -Z T i), (62) 
Jo i — 00 i % 1+1 

where P- lim denotes the limit in probability. This definition of stochastic integrals applies to the usual 
semimartingale setting where is a process in D and where Z is a semimartingalc. If Z = B H is a 
fractional Brownian motion process with Hurst coefficient H > |, the limit in (f6"2"| exists for a large 
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class of integrands, including continuous semimartingales and (^-functions of fractional Brownian mo- 
tion. In particular, the stochastic integral / B H dB H exists in the sense of (|62|) , and for the continuous 
semimartingale \P we have the following integration by parts formula, due to |39j : 

-J B H d^ + VB H = J ^dB H . (63) 

Before we state the main result of this section, we recall that a sequence {$ "} n gN of semimartingales 
defined on probability spaces (fi n , T n , P") is called good in the sense of if, for any sequence {Z n } ne ^ 
of cadlag adapted processes, the convergence >C-lim 71 _ >00 (^ r ™, Z n ) = (\f r , Z) implies the convergence 

C- lim ( V n ,Z n , [ Z n dV n ) = (v,Z, ( Zd$> 



Remark 4.1 Any semi-martingale ^> can be written as \& = M+A, in which M is a local martingale with 
Mq = and A is an adapted finite variation process. We will denote by [M, M] the quadratic variation 
of M and by \A\ the total variation of A. A sequence {^ n } n eN = M n + A n ) of semi-martingales is 
good on [0,T] if the sequences (E„{[M n , M n ]T}) n eN and (E n {\A\T}) n eN are bounded; see ]22l Theorem 
4-1] for details. 

In the following we will denote by {^"jngN a sequence of "good" semimartingales and by {Z n } n£ ^ a 
sequence of B- valued stochastic processes defined on some probability space (f2,.F, P) and assume that 
the following conditions are satisfied. 

Assumption 4.1 (i) The sample paths of the processes Z n are almost surely of zero quadratic vari- 
ation on compact sets, and ¥{Zq = 0} = 1. 

(ii) The stochastic integrals J^ n dZ n and J Z n dZ n exist in the sense of H62\), and the sample paths 
ti-> f*Z%_dZ? andt^ J*V2_dZ? are cadlag. 

We are now ready to state the main theorem of this section. Its proof requires some preparation which 
will be carried out below. 

Theorem 4.1 Let {$ "}„ £ n be a sequence of good semimartingales and let {Z n } n ^ be a sequence of 
Hi-valued stochastic processes that satisfy Assumption \4-l\ If ^ is a continuous semimartingale and 
if B H is a fractional Brownian motion process with Hurst parameter H > \, then the convergence 
£-lim n ^ 00 ('I' Tl , Z n ) — (?$l,B ) implies the convergence 



C- lim * 



Z n ,J ^ n dZ n ^j = (^>,B H , J ^dB H 



Before we turn to the proof of Theorem 14. li we consider an example where Assumption 14. II can indeed 
be verified. 

Example 4.1 Let {H n } ne ^ be a sequence of real numbers with H n > and assume that linin^oo H n = 
H > |. Let Z n be a fractional Brownian motion process with Hurst parameter H n and let ^ be a 
continuous semimartingale independent of Z n for all n. Since H n > i the processes Z n have zero 
quadratic variation. Moreover, C-liuin^oo Z n = B H because the centered Gaussian processes Z n and B H 
are uniquely determined by their covariation functions and all stochastic integrals exists in the sense of 
Thus, Theorem \4-l\ yields 

C- lim ( Z n , [ VdzA = ( B H , I ^dB H 



We prepare the proof of Theorem 14.11 with the following simple lemma. 
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Lemma 4.1 Under the assumptions of Theorem \4-l\ the processes [Z n ] and [Z n ,'9 n ] defined by 

[Z n ] t : = (Z t n ) 2 - 2 f Z n s _dZ n s and [Z n , * n ] t := - / Z n s _d^ n s - f V?_dZ?, 

Jo Jo Jo 

have P-a.s. sample paths which are equal to zero. 

PROOF. It follow from the representation of the stochastic integrals J Z n dZ n , J ^ n dZ n and / Z n d y S> n 
as probabilistic limits of Stieltjes-type sums that, for any t and each sequence of partitions {r'}; S N of 
[0, t] with lim/^oo max, |r- +1 - t\ |, 

[Z n ] = P-lim^(r, +i - Z^) 2 and [Z",*"] = P- hm ^ Z rf +1 ~ Z rO(Kl +1 - K0- 

i i 

Since a typical sample path of the stochastic integrals J ^> n dZ n and J Z n d^ n is in D, we can ap- 
ply the same arguments as in the proof of Theorem II. 6. 25 in |48j in order to obtain the inequality 
¥{[Z n , tf n ]? < [Z n ] t [*"] t } = 1. Thus, our assertion follows from F{[Z n ] t = for all t>0} = l. □ 

For the proof of Theorem 14. II we will also need the following result. 

Lemma 4.2 (i) Let C be the space of all real valued continuous functions. For n G N, let a n ,(3 n G D 
and assume that the sequence {(a n , /3 n )}neN converges in the Skorohod topology to (a, (3) £ CxC. 
Then, on compact intervals, the process 

In = (7n(*))t>o defined by j n (t) = a„{t)/3 n (t) 

converges to a/3 = (a(t)[3(t)) t >o in the Skorohod topology on ED. 

(ii) Let {(Y n , Z")}„ 6 n be a sequence of Hi-valued random variables defined on some probability space 
(0,^",P) that converges in law to (Y, Z). IfF{(Y, Z) G C x C} = I, then 

C- lim {{Y t n Zl l ) a < t < T } = (Y t Z t ) 0<t<T 

n — *oo 

holds for all T < oo. 

Proof. Since a and (3 are continuous, (i) follows from Lemma 2.1 in [35]. The second assertion 
follows from (i) and Skorohod's representation theorem. □ 

We are now ready to finish the proof of Theorem 14. II 



Proof of Theorem 14.11 Since {^™}„ eN is a sequence of good semimartingales and because a typical 
sample path of a fractional Brownian motion process is continuous, we deduce from Theorem 2.2 in [38] 
and from Lemma 14.21 (ii) that 



C- lim U> n ,Z n , J Z n dt> n ^j = (^£,B H ,J B H 'd&J and C- lim ($"Z n ) = {^>B H ) , 



(64) 



respectively. By the continuous mapping theorem, it follows from (|64p , from Lemma 14.11 and from the 
integration by parts formula for fractional Brownian motion ([53"]) that the finite dimensional distributions 
of the processes 



$> n ,Z n ,J y n dZ n j = \ V n ,Z n ,- J Z n dV n + f n Z r 
converge weakly to the finite dimensional distributions of the process 



V,B H ,- J B H dt> + VB H ^j = (^>,B H ,J ^dB H 



It also follows from that the sequence {/ Z"<Jf"} N is C-tight and the sequence {^ n Z n } n£N is 
tight. By Corollary VI. 3. 33 in J33j the sum of a tight sequence of stochastic processes with a sequence 
of C-tight processes is tight. Thus, continuity of the processes ^B H and J B H yields tightness of the 
sequence {/ W l dZ n } neW This shows that 



c- lim ( / ^ n dZ n I = ( B H , / y>dB H 
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□ 

In view of Lemma 14.11 and the integration by parts formula for fractional Brownian motion processes 
we also have the following approximation result for the integral of fractional Brownian motion process 
with respect to itself. 

Proposition 4.1 Under the assumption of Theorem \4-l\ it holds that 

C- lim / z:_dZ:= [ B?dB?. 
Jo Jo 



Proof. By Lemma |4~T1 



(Zl 1 ) 2 = 2 ( Z^_dZ 7 s l P-a.s. 
Jo 



Thus, in view of Lemma 14.21 (ii) and the Ito formula for fractional Brownian motion the sequence 
{(Z™) 2 }„ eN converges in distribution to 

ft 



(Bf) 2 = 2 / BfdBf. 
Jo 



This yields the assertion. □ 
We finish this section with the proof of Theorem 12.11 



PROOF of Theorem 12. 11 In terms of the processes Y and X £ introduced in ((Til) an( i UM> respectively, 
we have 



C- lim — ^=(X t e,Ar ) <t<T = C- lim 

JV^oo JN N^ao 



0<t<T 



V s Y s/e ds 

/ 0<t<T 

t 

*sdX E s 

/ 0<t<T 



and 



XUXt = I X e s Y s/£ ds 



10 jo 

where all the stochastic integrals have continuous sample paths. By Theorem 13. 1[ 

C- lim 1 y £ = cB H 

40 ei-»^L(e-i) 

with H > i and the sequence of semi-martingales consisting of a single element ^ is good as a result of 
Remark 14.11 since by assumption £'{[ V I', ^]t} < oo and E{|A|t} < oo. Therefore, the assertion follows 
from Theorem 14. II if we can show that the processes X e have zero quadratic variation on compact time 
intervals and that the stochastic integrals J ^dX e and J X e dX e exist as the probabilistic limits of 
Sticltjes-type sums. These properties, however, follow from (|2ip by direct computation. □ 



5. Conclusion & Discussion We proved a functional central limit theorem for stationary semi- 
Markov processes and an approximation result for stochastic integrals of fractional Brownian motion. 
Our motivation was to provide a mathematical framework for analyzing microstructure models where 
stock prices are driven by the demand of many small investors. We identified inertia as a possible source 
of long range dependencies in financial time series and showed how investor inertia can lead to stock price 
models driven by a fractional Brownian motion. Several avenues are open for further research. 

In our financial market model, investors do not react to changes in the stock price. Over short time 
periods, such an assumption might be justified for small, non-professional investors, but incorporating 
feedback effects, whereby traders' investment decisions can be influenced by asset prices, would be an 
important next step from an economic point of view. This, however, leads to a significant increase in the 
complexity of the dynamics, as discussed in Remark [ 
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Financial market models where stock prices are driven by the demand of many heterogenous agents 
constitute another mathematical challenge. Heterogeneity among traders has been identified as a key 
component affecting the dynamics of stock prices, and in recent years an extensive literature on the 
mathematical modelling of heterogeneity and interactions in financial markets has appeared. But the 
analysis is usually confined to models of active market participants trading every period, and it is therefore 
of interest to introduce heterogeneity and interaction into our model. One approach would be to build 
on the case study in Section [2.31 to look at models where inert traders (who alone would lead to a limit 
fractional Brownian motion price process) interact with active traders (who alone would lead to a limit 
standard Brownian motion price process) , and to include mechanisms by which their interaction leads to 
efficient markets. One could also try incorporate strategically interacting institutional investors; see [1] 
for a possible game theoretic framework. 

For long-term economic models, it is important to look at time- varying rates of long-range dependence, 
where the variation is caused by global economic factors, or regime changes. At some times, the market 
may be efficient, for example in a bullish exuberant economy like the late 1990s (see Figured]), while 
at other times, the Joseph effect may be prominent, such as during a recession or period of economic 
nervousness as in the early 1990s. The mathematical challenge is then to adapt the limit theorem of this 
paper to stochastically varying measures of inertia. 

Appendix A. The key renewal theorem in the heavy tailed case In this appendix we recall 
a result of Heath et al. ([E]) on the rate of convergence in the key renewal theorem in the heavy tailed 
case. 

Theorem A.l Let F be a distribution with domain [0, oo) satisfying 

F(t) = 1 - F(t) ~ t- a L(t) 

for some 1 < a < 2, and where L is a slowly varying function at infinity. Assume that F n is nonsingular 
for some n > 1. Let k = J Q F(x)dx be the expected value and denote by U the renewal function associated 
with F , that is, 

CO 

u = J2 pn - 

71=0 

Let z be a continuous, non-negative function of bounded variation on [0, oo), such that lim^oo z(t) = 0. 
That is, z(t) = J. C(dy) for some finite signed measure C on [0, oo). Let z* denote the total variation 
function of Q. That is, 

/CO 
\C\(dy). 

We will also assume that z is a directly Riemann integrable function (see p. 295 for a definition) on 
[0, oo), such that z(t) = o(F(t))and that 

z*{t) = 0(t- a+1 L{t)y (65) 

Let A = z(t)dt < oo. Then the function h : M + — * R + defined by 

A 

h(t) = -- z(t- s)U(ds) 
« Jo 

satisfies 

h(t) ~ - A t- a+1 L(t). 
[a — 1)k z 

Proof. The proof follows from the Remark on page 11 of [31]. □ 

Appendix B. Directly Riemann Integrable Functions The proof of our approximation result 
uses the notion of directly Riemann integrability. A real valued function g on R + is called directly 
Riemann integrable if 

CO CO 

lim y a inf g(t) = lim N a sup g(t) = 

These sums would be Riemann sums if not for the infinite limit of summation. If a function is directly 
Riemann integrable, it is Lebesgue integrable, but the converse is not necessarily true. 
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Lemma B.l (JTb\ Chapter 9], Proposition 2.16, (c) and (d)) 

(i) Let g > be a monotone non-increasing junction; then g is directly Riemann integrable if and 
only if g is Riemann integrable. 

(ii) Let g > and let <b be a distribution function on R + . If g is directly Riemann integrable, then 
4> * g is directly Riemann integrable as well. 
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